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Ti^ij- ' Abstract 

I Weakly-interacting massive particles (WIMPs) annihilating in the center of the Sun or the Earth 

. may give rise to energetic neutrinos which might be discovered by astronomical neutrino detectors. 

' The angular distribution of the neutrino-induced muons is considered in detail via Monte Carlo 

\ simulations. It is shown that large underground Cerenkov neutrino telescopes might be able to 

extract the WIMP mass from the width of the muon angular distribution. 

Oh! 

^ ■ 1 Introduction 



Weakly-interacting massive particles (WIMPs) with masses in the GeV-TeV range are among the leading 
non-baryonic candidates for the dark matter in our galactic halo. One of the most promising methods for 
the discovery of WIMPs in the halo is via observation of energetic neutrinos from annihilation of WIMPs 
in the Sun and/or the Earth Q. Through elastic scattering with the atomic nuclei in the Sun or the 
Earth, a WIMP from the halo can lose enough energy to remain gravitationally trapped Trapped 
WIMPs sink to the core of the Sun or the Earth where they annihilate into ordinary particles: Icptons, 
quarks, gluons and - depending on the masses - Higgs and gauge bosons. Because of absorption in the 
solar or terrestrial medium, only neutrinos are capable of escaping to the surface. Most WIMP candidates 
- among them the supersymmetric candidate, the neutralino - do not annihilate into neutrinos directly 
Nevertheless energetic neutrinos are eventually produced via hadronization and/or decay of the 
annihilation products. These energetic neutrinos may be discovered by astronomical neutrino detectors. 

Trapped WIMPs are (to a good approximation) in thermal equilibrium with the core of the Sun 
and/ or the Earth. The radial extension of the WIMP annihilation region is a function of the WIMP mass 

||], heavier WIMPs lying deeper in the core. This led Gould (?) to suggest that the WIMP mass might 
be inferred from the angular size of the annihilation region. 

In this letter, we consider Cerenkov neutrino telescopes. They consist of large underground arrays of 
photo-multipliers to detect the Cerenkov light emitted by muons generated in charged-current interactions 
of neutrinos with the medium surrounding the detector. Underground Cerenkov detectors, originally built 
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to search for proton decay, have aheady started to explore (and constrain) WIMP dark matter candidates 

Here we want to study if Cerenkov neutrino telescopes currently planned or under construction Q 
might realistically expect to be able to extract the WIMP mass from the muon angular distribution 
(once measured). We include the uncertainties in the determination of the neutrino direction due to 
the neutrino-muon scattering angle in charged-current interactions, to multiple Coulomb scattering of 
the muon on its way to the detector and to an intrinsic angular resolution in the determination of the 
direction of the muon track. 



2 Annihilation channels and muon fluxes 

WIMPs trapped in the core of the Sun and/or Earth can annihilate to a fermion-antifermion pair, to gauge 
bosons, Higgs bosons and gluons (xx l+l^, qq, gg, qqg, W+W-, Z^Z^, Z°H°, W^H^, H°H°). These 
annihilation products will hadronize and/or decay, eventually producing high energy muon neutrinos. 

Edsjo reconsidered the whole chain of processes from the annihilation products in the core of the 
Sun or the Earth to detectable muons at the surface of the Earth. He performed a full Monte Carlo 
simulation of the hadronization and decay of the annihilation products using Jetset 7.3 of the 
neutrino interactions on their way out of the Sun and of the charged-current neutrino interactions near 
the detector using Pythia 5.6 and finally of the multiple Coulomb scattering of the muon on its 
way to the detector using distributions from Re f. | l2 | . 

With respect to previous calculations jl^, |l4|, the Edsjo Monte Carlo treatment of the neutrino 
propagation through the Sun bypasses simplifying assumptions previously made, namely neutral currents 
are no more assumed to be much weaker than charged currents and energy loss is no more considered 
continuous. In the new treatment, the neutrino energy spectrum at the surface of the Sun is obtained as 
follows. The thickness of the Sun and the neutrino mean free path determine the probability of neutrino- 
nucleus interactions. Each interaction is randomly chosen to be a charged-current interaction, in which 
case the neutrino is considered absorbed, or a neutral-current interaction, in which case the neutrino 
is degraded in energy according to distributions in Pythia 5.6. The procedure is continued until the 
neutrino has reached the surface of the Sun. The resulting neutrino spectrum differs significantly from 
previous calculations only in the high energy tail. But from this high energy tail comes the most important 
contribution to the muon flux in Cerenkov neutrino detectors.^ Hence at a WIMP mass of 1500 GeV (50 
GeV) Edsjo finds a muon flux 20% (5%) higher than that obtained by Ritz and Seckel fl^ . 

For more details on Edsjo results, we refer the reader to Ref. In the following, we rely on his 

results obtained by simulating 10^ WIMP annihilation events per annihilation channel and WIMP mass. 



3 Annihilation profiles 

The annihilation region in the Sun can be regarded as point-like, its angular size being < 0.005° for the 
WIMP masses we are interested in, m > 10 GeV. For the Earth, on the contrary, the angular extension of 
the annihilation region is non- negligible and decreases as one over the square root of the WIMP mass Q . 
In fact, the annihilation rate per unit volume at a distance r from the center of the Earth is proportional 
to the square of the WIMP number density n{r). The latter may be written as 

n(r) = n(0)e-'''/2'-x, (1) 

with 

3kT y^\_ 0.56i?e (2) 
ATTGpm] ^m/GeV' 

Here we have taken the radius of the Earth i?^ ~ 6400 km, the central Earth temperature T ~ 6000 K 
and the central Earth density p ~ 13 g cm^'^. For the WIMP masses we are interested in, it is a very 
good approximation to consider a constant Earth density in the region where WIMPs are concentrated. 

^We remind that this is so because the muon flux is the product of the neutrino flux by the charged-current cross section 
and the muon range, and both of these are proportional to the neutrino energy. 
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Figure 1: Projected angular distributions of WIMP-generated neutrinos from the Earth for WIMP masses 
of (solid line) 50 GeV, (dashed line) 100 GeV, (dotted line) 200 GeV, (dash-dotted line) 350 GeV and 
(wide dotted line) 750 GeV. The analogous distributions from the Sun are simply a narrow peak atOu — O. 



For an observer close to the surface of the Earth the angular distribution of the neutrinos generated in 
WIMP annihilations results 

H <E> 2 2 

, , " (xe-° '^'^ erf (a cos 0.) (3) 
d<y9^dcos 

where is the angle between the neutrino direction and the center of the Earth, is the associated 
azimuthal angle, a — R^/r^ and erf is the error function. 

This expression simplifies in our case, m > 10 GeV, for which 6^ is typically smaller than 15°. In 
this case the 0i,-distribution is simply approximated as a gaussian in the transverse plane, the plane 
orthogonal to the directions of either the Sun or the Earth centers, 

oce-^(«^+««), (4) 



d^.d^j, 



where Ox and Oy have obvious meaning. 

To reduce the fluctuations due to limited statistics it is more convenient to consider the projected 
distribution in Ox- From Eq. (^) we find that the projected distribution can be considered gaussian. 



d£, 
AOx 

with root mean square value 



(5) 



1 r 23° 
^ ^rad ^— — (m > 10 GeV) . (6) 
V2 ^TO/GeV 

Projected angular distributions of the neutrino flux from the Earth arc shown in Fig. |l| for WIMP 
masses between 50 and 750 GeV. The distribution width decreases for increasing WIMP masses. Note 
that these distributions are independent of the neutrino energy spectrum and of the specific annihilation 
channel. The analogous distributions for the Sun are simply narrow peaks at On = 0. 



4 Muon angular distributions 

In Cerenkov neutrino telescopes it is not possible to measure the angular distribution of the neutrinos 
directly since it is the muon produced in charged-current interactions that can be detected. The direction 
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Figure 2: Projected angular distributions of neutrino-induced muons from WIMP annihilations in (a) the 
Earth and (h) the Sun for WIMP masses of (solid line) 50 GeV, (dashed line) 100 GeV, (dotted line) 200 
GeV, (dash-dotted line) 350 GeV and (wide dotted line) 150 GeV. The distributions are shown for hard 
channels (W~^W~ for 100-750 GeV and t^t^ for 50 GeV), with a detector muon threshold E*^ = 10 
GeV and a detector angular resolution Odet — 1-4°. 



of the neutrino is somewhat lost because of two effects: (1) the muon produced in a charged-current 
interaction exits at an angle Oqc with respect to the incoming neutrino and (2) the same muon undergoes 
multiple Coulomb scattering on its way to the detector, changing direction by an angle Ocoui- Both angles 
are approximately gaussian in the transverse plane (at least in the central region, for 9cc has non-gaussian 
tails) with root mean square values 



nrms ^ 
^CC 



y/EjGeV 
and 



nrms 
^Coul 



(E^ > 10 GeV) (7) 



O 1 o 

(8) 



y%jGeV 

where the first relation is obtained by simulations with Pythia 5.6 and the second relation is from 
Ref. Notice that both angles get smaller with increasing neutrino (and muon) energy. 

There is an additional uncertainty coming from the reconstruction of the muon track. Each neutrino 
telescope has an intrinsic angular resolution in determining the direction of the muon. We assume that 
the error in its determination is normally distributed with root mean square value Odet, typically of the 
order of 1°. 

In Fig. ^ we plot the muon angular distributions for hard channels in the Earth and in the Sun, 
obtained from the full Edsjo simulations in |l^. These distributions are representative of any hard 
neutrino spectrum. For a neutrino spectrum to be hard it is not necessary that it is dominated by a hard 
channel, like W'^W~, Z'^Z^ and t~^t~. Because of the previously-mentioned importance of the high 
energy tails, it suffices that the branching ratio into hard channels is greater than ^ 10%. Softer neutrino 
spectra, e.g. those from the bb and H^H'^ channels, give rise to wider angular distributions. This is due 
to the energy dependence of 6'^™'' and fi'p'oui ■ Note the difference between the neutrino (Fig. |]) and muon 
(Fig. H) angular distributions: charged-current interactions and multiple Coulomb scattering make the 
width dependence on WIMP mass stronger. 
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Figure 3: WIMP mass m versus full width half maximum 0fwhm of the neutrino-induced muon distri- 
bution for soft (hh) and hard (W'^W~ and t'^t~) channels. The solid lines correspond to Ej^ = 10 GeV 
and the dashed lines to Ej^ = 2 GeV. a) For Sun-bound WIMPs with a detector angular resolution 
of Odet =^ 1.4°. b) For Earth-bound WIMPs with a detector angular resolution of 9det = 1.4°. c) For 
Earth-bound WIMPs with a perfect angular resolution. 

5 WIMP mass determination 

Information on the WIMP mass might be extracted from the width of the distribution in the transverse 
plane. With Hmited statistics it might be convenient to project this two-dimensional distribution onto 
the 9x axis. As a measure of the width we choose the full width half maximum OpwHM of the projected 
distribution. The FWHM is not sensitive to the non-gaussian tails of the distribution, which reflect the 
non-gaussian tails in dec- Moreover, as long as the detector resolution 9det is small with respect to 
6 FWHM, it should be relatively easy to extract the FWHM even in the presence of a background (muons 
and neutrino-induced muons from cosmic ray interactions in the atmosphere). 

In Figs, ^-c we show the dependence of the WIMP mass on the full width half maximum for some 
representative cases. We present the soft 55-channel and the hard W^W~- and r+r^-channcls for the 
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Figure 4: An example of mass determination from a projected muon angular distribution. The simulated 
histogram includes the expected background in one year of exposure with a 1 km? detector and a typical 
signal from annihilation of 100 GeV WIMPs in the Earth. The solid line is a fit of a gaussian plus a 
constant. 



Earth and the Sun for two different muon energy thresholds, = 2 GeV and Ejl" = 10 GeV. Figs. §i 
and ||b include a detector angular resolution 9det = 1-4°. For the sake of comparison, Fig. |3|3 shows 
the ideal case of a perfect angular resolution. We have checked that all distributions are indeed well 
approximated by gaussians in the central regions. 

We see that the detector angular resolution is the limiting factor for the mass determination of heavy 
WIMPs (m > 400 GeV). For lighter WIMPs, it seems promising to infer their mass from the Earth 
muon distributions. We remind that this is also the mass range in which the signal from the Earth is 
expected to be significant (|] . The WIMP mass could also be extracted from the Sun muon distributions 
provided the detector muon energy threshold is low. In fact, the width of the angular distribution for the 
Sun is dominated by charged-current scattering, which acts as a mass spectrometer in diffusing neutrinos 
according to their energies and so according to the WIMP mass. 

From Figs. and analogous ones, approximate relations between the error on the WIMP mass 

and the error in Ofwhm may be obtained. For example, for E'^'' — 10 GeV, 9det = 1-4° and m < 300 
GeV, we read from Fig. |^ 

Alogm- 0.15 + 0.15A6ii.vi/HM (9) 

with /S.9pwHM in degrees. The first term represents our ignorance on the actual annihilation channel 
and the second term comes from the slope of the m - OpwHM curve. For AOpwHM ~ 1°, relation 
gives an uncertainty on the WIMP mass of roughly a factor of 2.0. 

We present now an example of mass determination for Earth-bound WIMPs. We consider one year of 
exposure of a 1 km^ detector with a muon energy threshold of 10 GeV. The atmospheric background in 



the direction of the center of the Earth is expected to be 20 muons per square degree |15|. We choose a 
WIMP mass of 100 GeV and generate a signal of 2000 muons, a reasonable number for supcrsymmetric 
models with neutralinos of this mass. In this way we obtain a simulated muon angular distribution in a 
15° X 15° region centered towards the center of the Earth. We then analyze these simulated data. We 
project the muon distribution onto the Ox axis and obtain the histogram shown in Fig. H By fitting a 
gaussian plus a constant to this histogram we obtain a full width half maximum Ofwhm — 8.9° ± 0.7°. 
From Fig. we then read the WIMP mass range corresponding to the fitted FWHM range, m — 90^25 
GeV. We are satisfied that the mass range obtained contains the original WIMP mass. The uncertainty 
is approximately a factor of 1.5, in agreement with the approximate relation Eq. (^. One might worry 
about the size of the uncertainty, but we believe that even a rough determination of the WIMP mass 
would be of enormous importance. We are exploring ways to reduce this uncertainty. 
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In favorable cases, many muon neutrinos might be detected from both the Earth and the Sun. Hav- 
ing both angular distributions, the angular smearing due to charged-current interactions and multiple 
Coulomb scattering would be directly represented by the Sun muon distribution (recall that the annihila- 
tion region in the Sun can be considered pointlike). This distribution might then be subtracted from the 
Earth muon distribution, leaving the neutrino angular distribution from the Earth. From it the WIMP 
mass could be obtained in a direct way via Eq. (||). Notice that for heavy WIMPs (to > 100 GeV) one 
should correct for absorption of neutrinos on their way out of the Sun, e.g. by specifying the shape of 
the neutrino energy spectrum. 

6 Conclusions 

WIMPs annihilating in the center of the Sun or the Earth may give rise to a neutrino-induced muon 
flux in astronomical neutrino detectors. The width of the muon angular distribution carries information 
on the WIMP mass, because the size of the annihilation region, the charged-current neutrino- nucleon 
scattering and the multiple Coulomb scattering of the muons all depend on the WIMP mass. Detailed 
Monte Carlo simulations have been used to obtain the muon angular distribution for WIMP annihilations 
in the Earth and in the Sun. It has been shown that the WIMP mass can be inferred, for WIMPs lighter 
than ^ 400 GeV, from the Earth distribution and, provided the muon energy threshold is low (< 5 GeV), 
also from the Sun distribution. This seems therefore a promising method of determining the WIMP mass 
and we look forward for the detection of a WIMP signal in neutrino telescopes. 
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